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Recently, properties of collective states of interacting non-abelian anyons have attracted a con- 
siderable attention. We study an extension of the 'golden chain model', where two- and three-body 
interactions are competing. Upon fine-tuning the interaction, the model is integrable. This provides 
an additional integrable point of the model, on top of the integrable point, when the three-body 
interaction is absent. To solve the model, we construct a new, integrable height model, in the spirit 
of the restricted solid-on-solid model solved by Andrews, Baxter and Forrester. The heights in our 
model live on both the sites and links of the square lattice. The model is solved by means of the 
corner transfer matrix method. We find a connection between local height probabilities and char- 
acters of a conformal field theory governing the critical properties at the integrable point. In the 
antiferromagnetic regime, the criticality is described by the Zk parafermion conformal field theory, 
while the ''"(^)i'"'"(^)i''^"(^)fc-2 goset conformal field theory describes the ferromagnetic regime. 
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I. INTRODUCTION 

The last half decade has seen a big increase in the interest of topological phases of matter. In this paper, we will 
study a model which is inspired by the prototype of a topological phase, namely the (fractional) quantum Hall effect. 
It has been conjectured that there exist fractional quantum Hall states with excitations which exhibit non-abelian 
statistics [T]. One of the key properties of this type of excitations, called non-abelian anyons, is that a topological 
state with a number of non-abelian anyons present, is degenerate. The number of degenerate states is exponential, 
while the energy splitting in real systems decays exponentially with the average distance between the non-abelian 
anyons. 

An important question which raises itself is what happens if the anyons are close to one another, such that they 
start interacting. To this end, a one-dimensional (ID) model of interacting anyons was constructed in Ref. [2 , called 
the 'golden chain', because it was based on Fibonacci anyons. The philosophy behind this model was to stay as closely 
as possible to a Heisenberg model of interacting spins. In fact, the golden chain is precisely that, a Heisenberg model 
with two-body nearest neighbor interactions, but for anyons instead of spins. The phase diagram of such anyonic 
Heisenberg models turns out to be rich, even richer than the phase diagrams of the ordinary spin case. We will not 
embark on a long discussion of the phase diagrams of these models here, but focus in the next section on one particular 
example of interest for the current paper, an extension of the golden chain model with competing two- and three-body 
interactions [3]. We would like to point out that studying the effects of interacting anyons in ID gives insight into 
the fate of interacting anyons in two-dimensional (2D) systems. The interactions between the anyons can nucleate a 
new topological liquid, and the collective behavior of the ID chain describes the boundary between the original and 
nucleated topological phases, see Refs. [HIS] and Refs. [SJ[7] for related work. 

An interesting property of the anyonic chain models is that they exhibit (fine-tuned) points, at which it is possible to 
solve the model exactly. Obviously, having access to an exact solution, even though such a solution is only available at 
special points, greatly enhances the understanding of the model. The golden chain (with a two-body interaction), and 
its cousins which are obtained by replacing the Fibonacci anyons by anyons based on su(2)fc, can be mapped onto a 
two-dimensional classical statistical mechanics model, namely the restricted solid-on-solid (RSOS) model, introduced 
and solved by Andrews, Baxter and Forrester (ABF) 0. This RSOS model consists of heights living on the sites of 
the square lattice. Plaquettes are weighted depending on the heights of the sites forming the plaquettes. To solve the 
model, ABF employed the so-called corner transfer matrix (CTM) method It was found that the model exhibits 
various ordered phases, separated by critical points. The connection between the critical exponents and conformal 
field theory (CFT) was made by Huse [TD]. The anyonic chains correspond to the RSOS model at the critical point, 
and are therefor critical themselves, and governed by the same CFT [2,. 
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The CTM method allows one to calculate, in the limit of infinite lattice size, the probability for a site in the 
bulk to have a particular height. Interestingly, it has been observed that off-critical local height probabilities of 
integrable models are intimately related to partition functions of the associated critical theories in a finite box with 
appropriate boundary conditions Moreover, close to critical points, these height probabilities are given in terms 
of characters of a CFT, which describes the critical behavior of the model [12J. The characters stemming from the 
height probabilities in the RSOS model, are for instance given in Ref. [13]. These characters can be interpreted in 
terms of fractional exclusion statistics [T^ , or more specifically, a non-abelian version thereof |151 In light of the 
current paper, we would like to point out that the height probabilities one obtains for a finite system, correspond to 
finitized characters. Interestingly, the opposite 'ends' of these finitized characters correspond to the full characters of 
two different conformal field theories, in the large size limit. These two different conformal field theories describe the 
critical properties of different critical points of the (RSOS) model. In terms of the anyonic chains, these two critical 
points are related to each other by changing the overall sign of the interaction. 

In this paper, we will examine an integrable point of anyonic chains with competing two- and three-body interaction 
terms. To solve the model, we introduce a new statistical mechanics model, which builds on the RSOS model, which 
has six different types of plaquettes. In our new model, we will combine four plaquettes of the RSOS model, and shift 
two of the plaquette weights, to obtain a non-trivial generalization. Via this procedure, one obtains a model with 
sixty-six different types of plaquettes. This composite height model gives rise to the anyonic chain Hamiltonians of 
interest, via the usual 'anisotropic limit'. To solve the model, we will follow the work of ABF, and employ the corner 
transfer matrix method to calculate the height probabilities. 

A very closely related loop model has been studied in the literature |17l \TE\ . In fact, the i?-matrix, used to 
construct the row-to-row transfer matrix, has the same underlying algebraic structure, namely that of the Temperley- 
Lieb algebra, and the construction of our 2D model is inspired by the work in Refs. |17l I18| . However, the quantum 
chain of that work is defined on a completely different Hilbert space, in comparison to the anyonic chains. In our 
case, the Hilbert space does not have a tensor product decomposition, which hinders solving the model by means of 
the Bethe Ansatz, which is the method used in Refs. |171 [18] (in their representation, the Hilbert space does have a 
tensor product decomposition). The Hilbert spaces of the anyonic chains exhibit a non-local, topological symmetry, 
giving rise to a topological quantum number, which can be used to label the eigenstates of Hamiltonians respecting 
this symmetry (see, for instance Refs. [21 En])- In gapped phases where the ground state breaks this symmetry, one 
finds additional degeneracies. This additional structure, which is intimately tied to the non-abelian nature of the 
anyons, seems to be absent in the work of Refs. [TT] [T5]. 

The outline of the paper is as follows. In Section |TT] we briefiy introduce the anyonic chain Hamiltonians, starting 
with the original golden chain, and its generalization by introducing the three-body interaction. We also briefiy discuss 



the anyonic chains based on su{2)k anyons. In Section III we discuss the connection of the golden chain with the 



RSOS model. The following Section [IVj contains the definition of our new, composite height model, which builds on 
the RSOS model. The corner transfer matrix method is described in Section |V] which leads to the expressions for the 
height probabilities. These are used in the Section VI to obtain information about the various phases of the model. 
In Section |VIH we calculate off-critical local height probabilities in different phases and show that they are given in 
terms of characters of the CFTs, which govern the corresponding critical points. We conclude in Section |VIII[ In 
Appendix |X] we explicitly give the different types of plaquettes of the composite height model. Appendix [B] contains 
various limits of the CTM's, which are used in the main text. Appendix [C] deals with a certain limit of the plaquettes. 
Finally, in Appendix |D] we give the details of the connection of the various heigh probabilities with CFT characters. 



II. ANYON CHAINS WITH COMPETING INTERACTIONS 



A. The golden chain 

The first model of interacting anyons, introduced in Ref. [2], was dubbed the 'golden chain'. In this model, so- 
called Fibonacci anyons interact in basically the same way as spins in the Heisenberg spin chain, namely, an energy 
is assigned depending on the overall spin state of two interacting particles. In the S = 1/2 Heisenberg chain, the 5* = 
state of two neighboring spins is favored for antiferromagnetic interactions, while for ferromagnetic interactions, the 
5 = 1 channel for neighboring spins is favored. 

To explain the interaction of the golden chain, we will first briefiy introduce the notion of (Fibonacci) anyons. A 
rather extensive introduction on this topic can be found in Ref. [W. The Fibonacci anyon model consists of two 
types of particles: the 'trivial' or vacuum particle 1, and the Fibonacci particle r. As with ordinary spins, one can 
combine, or fuse these particles, and decompose the product. This fusion product is the direct analog of taking tensor 
products of spins. Contrary to su{2) spin, there is only a finite number of types of particles. In addition, there is no 
internal quantum number, such as s^. The reason for this will become clear shortly. 
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The rules for combining the anyons in the Fibonacci model are as follows 1x1 = 1, Ixt = t and t x r = 1 + r, the 
latter being the only non-trivial fusion rule. 

Let us take a chain of L t anyons. To describe the Hilbert space of this system, it is easiest to think in terms of a 
so-called fusion chain, as depicted in figure [T] 

T T T T T 



Xo Xi X2 X2, 

Figure 1: The fusion chain consisting of Fibonacci anyons. Consistent labelings {xo,xi, ...,xi^) form the Hilbert space of the 
anyonic chain. 

This fusion chain consists of labelled lines, L incoming lines labelled r, which represent the Fibonacci anyons that 
form the chain. The lines connecting these incoming r anyons are labelled xq, xi, etc. These labels, which can take 
the values 1 and r, are the 'degrees of freedom'. The set of consistent labelings {xo,Xi, . . . ,xl) forms the Hilbert 
space of the chain. For a labeling to be consistent, the fusion rules have to be satisfied at every vertex. This means 
that Xi+i has to be in the fusion of Xi x r. This means that one can not have Xi = Xi+i = 1, because this would 
violate the fusion rule 1 x r = r. Apart from the constraint that no two neighboring labels can both take the value 
1, the labelings are arbitrary. Because of the constraint, the size of the Hilbert space grows as , where d < 2. 
It is not so hard to convince oneself that in fact d is the golden ratio, d = f = {1 + \/5)/2. In the remainder of the 
description, we will assume periodic boundary conditions, xq = xl- In this case, the size of the Hilbert space is given by 
dimHi = Fib(L+ 1) +Fib(L- 1), where Fib(n) is the n^^ Fibonacci number, defined by Fib(n) = Fib(n- 1) +Fib(n- 2) 
and the initial conditions Fib(O) = and Fib(l) = 1. Loosely speaking, one can say that each Fibonacci anyon has 
a fractional number of degrees of freedom, namely d, explaining the absence of an internal quantum number. More 
importantly, one can not assign a local Hilbert space to each anyon. This is the reason we had to resort to the fusion 
chain to describe the Hilbert space, which can not be described as a tensor product of local Hilbert spaces, as is the 
case for ordinary spin chains. 

We would like to point out that in the description of the Hilbert space, we did not make use of the braid properties 
of the Fibonacci anyons. Often, the braid properties are used to define the concept of non-abelian statistics. What 
we have done here instead, is to use the fusion properties of the Fibonacci anyons, which make non-abelian statistics 
possible. In particular, in order for non-abelian statistics to be possible, one needs a Hilbert space whose dimension 
is at least two. This in turn is possible, if one considers particles (anyons) which have multiple fusion channels upon 
fusion with another particle, such as r x t = 1 + r. In the construction of interacting anyonic chains, it is the presence 
of multiple fusion channels which is the key property of non-abelian statistics which is utilized. One can consider 
models in which the explicit, non-abelian braid properties are used to define the Hamiltonian (see for instance [2(J|), 
but that is not the route will take in the present paper. 

We turn our attention to the description of the Hamiltonian. We first concentrate on the Hamiltonian of the original 
golden chain model. The interaction between two anyons depends on their overall fusion channel. Favoring the overall 
fusion channel of two neighboring anyons to be 1 will be called antiferromagnetic, while favoring the r channel will 
be called ferromagnetic interaction. 

In our description of the Hilbert space in terms of a fusion chain, the fusion channel of two neighboring anyons is not 
explicit, because the lines associated with these two anyons do not meet in one vertex. One can, however, perform a 
local basis transformation, which makes this fusion channel explicit. The matrix describing this basis transformation 
is called the i^-matrix, which is the direct analog of the 6j-symbols in the case of su{2) spin. These describe the 
change of basis between the following two possible ways of describing the Hilbert space of three spins: {Si ® S2) ® 
and 5*1 <S> (5*2 <S> S'3). In figure |2] we depict the _F-matrix in terms of the fusion-chain pictures. 

Figure 2: The F-matrix elements as a local basis transformation. 

After performing this local basis transformation, the fusion channel of the two neighboring anyons is explicit and 
it is given by x'^. Based on the value of this label (which can be 1 or r), we can now assign the energy by projecting 
onto the 1 or t channel. All that is left to do is to perform one more basis transformation, to go back to the original 
basis. 
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To make the model completely explicit, we have to give the values of the F-matrix elements. For the anyon 
models based on su{2), there exist explicit formulas for the i^-matrices, see for instance Ref. [3T]. We will not go 
into the details here, but simply give the results for the F-symbols we need. Note that there is some degree of 
freedom in the i^-symbols, the explicit form we give fixes this. The F-symbols we need are specified by specifying 
the values of {xi-i,Xi,Xi+i), which can take the values {(1, r, 1), (1, r, t), (r, r, 1), (r, 1, r), (r, r, t)}. In addition, 
we have to specify the degrees of freedom after the basis transformation, {xi^i,x[,Xi+i), which can take the values 
{(1, 1, 1); (1, r, r); (r, T, 1); (r, 1, r), (r, T, r)}. Using this ordering of the states, the F-matrix takes the form (high- 
lighting the block structure) 



F 



n 



\ 



(1) 



The inverse transformation is given by the same matrix, because F = F ^ , as is easily checked by using (f = 1 + d. 
We can now easily form the local projection operators -Pj^^^^y and P2^l,Q^y, which project onto the r and 1 channels, 
i.e. give energy to these channels. Hence, F^^I^Q^jy corresponds to the antiferromagnetic interaction. Explicitly, 



these projection matrices take the form -Pgl^ody ^ ^ ' diag(0, 1, 1, 0,1) ■ F and ^2^^^^^ = F ■ diag(l, 0, 0, 1,0) • F. The 

components of these matrices read {F!1:^^^'''^''^Y_^'' [Fl^^-^''^''^y^ and {Fxi+i''^''^Yi {Fxi+i'''^''^)^. ' respectively. ExpHcitly 
written out, this becomes 



(1) 



2-body 



/I 



d-2 



\ 



d- 
d-'l 



r,(r) 
2-body 



/o 



-d-i 



\ 



-d-i 
^2 



(2) 



One can easily check that Fg^^^^y + F2l^ody ~ ^' where 1 is the identity matrix. 



{r) 



We can now write down the golden chain Hamiltonian as the sum of the projection operators ^'2r^odyi' 

J2 Zi=i ^2-body i' whcrc the projector F2_,^Qjy ^ assigns a positive energy if the anyons i and i + 1 are in the r channel. 
In the original golden chain paper |2 , it was shown numerically that this model is critical for either sign of the 
interaction. Moreover, the central charge was determined via the entanglement entropy, resulting in c = 7/10 and 
c = 4/5 for antiferromagnetic ( J2 = 1) and ferromagnetic ( J2 = -1) interactions, respectively. Exact diagonalization of 
the model showed that for antiferromagnetic interactions, the low-lying part of the spectrum can be described in terms 
of the minimal model A^4.5, describing the tri-critical Ising model. The critical model describing the ferromagnetic 
system is that of the Z3 parafermions. 

It was subsequently realized that the golden chain Hamiltonian can in fact be obtained from an exactly solvable 
model, a particular version of the restricted solid-on-solid models |8j. These models exhibit various ordered phases, 
separated by critical points. It is these critical points the golden chain can be mapped to. As a result, one can obtain 
information about the critical theory of the golden chain, by studying the critical behavior of the RSOS models. We 
will discuss these RSOS models, and their connection to the anyonic chain Hamiltonians in more detail in the next 
section. 



H 



B. Competing interactions 



Having introduced the golden chain Hamiltonian, in which the anyons interact via a two-body nearest-neighbor 
interaction, we now consider the effect of introducing a three-body interaction, which was first considered in Ref. [3]. 

It is well known that if one adds a three-body term (with large enough coupling) to the S = 1/2 Heisenberg 
antiferromagnet, a gap opens, and one enters the Majumdar-Ghosh (MG) phase [22, 23!. The phase diagram of the 
Heisenberg chain with competing nearest-neighbor two- and three-body interactions is rather rich. Similarly, it was 
expected that the phase diagram of the golden chain with competing two and three-body interactions is rich as well. 
This model was studied in Ref. |3j, indeed finding an interesting phase diagram (see Fig. |4]), which we briefly review 
below, after introducing the details of the three-body interaction. 



Figure 3: The F transformations needed for the three-body interaction. 



In order to find the fusion channel of three neighboring anyons, we first have to perform two F transformations, 
after which this fusion channel is explicit. One can then project onto the desired channel, and go back to the original 
basis. This was explained in detail in Ref. [20' . The schematics of the basis transformation is given in figure [s] 

The three-body interaction will depend on four labels {xi-i^Xi,Xi+i,Xi+2)- To give the interaction matrix, we use 
the following basis {(1, r, r, 1); (1, r, 1, r), (1, r, r, r); (r, 1, t, 1), (r, t, t, 1); (r, 1, r, r), (r, r, 1, r), (r, r, t, r)}. In this 
basis, the three-body projectors take the following form 



/I 



r,(l) 

3-body 














d-2 

-d-l 



-d-l 



-d-i 

-d-i 



/o 



3-body 



d-l 



d-l 



d-i 
d-l 
2d-2 



(3) 



In terms of these projectors, the most general interaction we can write down takes the form 



i=l 



lOP. 



(r) 



3-body, 2 ■ 



(4) 



In this equation, the three-body projectors act on the quadruples {xi-i,Xi,Xi+i,Xi+2), while the two-body projectors 
act on triples {xi-i,Xi,Xi+i). Moreover, we assume periodic boundary conditions, xi+l = Xi, and we introduced the 
couplings J2 = cos 9 and J3 = sin 6. 



C. Phase diagram of the J2 J3 model 



We now briefly describe the phase diagram of the J2 - J3 model, as a function of the angle 9. For more details, we 
refer to Refs. [31 [53]. The phase diagram of this model is shown in figure |4] 




Figure 4: Phase diagram of the J2 - J3 model; J2 = cosO, J3 = sin^. 
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The angle 9=0 corresponds to the original, two-body golden chain model, which is critical and described by the 
tri-critical Ising model. For a finite range of both positive and negative angles, this behavior persists, so we have an 
extended range of critical behavior. At the angle given by tan 6* = with both the two- and three-body interactions 
antiferromagnetic, there is a phase transition to a gapped phase. In this gapped phase, the ground state is two-fold 
degenerate, with all the ground states occurring at zero momentum. It turns out that at the transition point, which 
is described by the parafermion theory, a non-local symmetry, dubbed topological symmetry, is broken. 

The critical phase around 9 = also gives way for a gapped phase if 9 is decreased, namely around 9 « -0.4727r. In 
the resulting gapped phase, both the spatial and topological symmetries are broken, giving rise to four (dimerized) 
ground states. This phase is the anyonic equivalent of the Majumdar-Ghosh phase in spin-1/2 chains with a purely 
three-body interaction. The phase transition from the tri-critical Ising region to the MG-like phase is described by 
the tetra-critical Ising model. 

The critical behavior at 9 = n, the original golden chain with ferromagnetic interactions is described in terms of 
the Z3 parafermions (the critical behavior of the 3-state Potts model). This point is part of an extended critical 
region, which extends all the way to the gapped phase in the region when both couplings are antiferromagnetic. This 
endpoint of that gapped phase is around 9 « O.SIGtt. 

The other end point of the extended critical region containing the point 9 = n is marked by a first order transition 
located at « 1.0757r to a sliver of an incommensurate region, which quickly gives way to an extended critical region. 
This region has low-lying states at momenta K = 0,7: 12,7: . in 12, and was therefor dubbed the Z4 phase. This phase 
has a transition to the MG phase around 9 ^ I.SStt. 

We will close this quick walk through the phase diagram by noting that there are two special points which lie in the 
gapped phases, namely at tan0 = 1/5/2 and 9 = 3tt/2. At these special points, the ground states of the gapped phases 
are exactly degenerate (as opposed to exponentially degenerate with system size) , and moreover, one can obtain these 
ground states explicitly. For more details about the phase diagram and its peculiarities, we refer to Refs.[3, 24J. 

In this paper, we will mainly concentrate on a special angle, given by tan0 = 1/ip, which corresponds to (on the 
one hand), the transition between the extended critical region at antiferromagnetic two-body interactions and the 
gapped phase obtained by introducing the antiferromagnetic three-body interaction. Upon changing the sign of the 
Hamiltonian, one ends up in the extended Z4 critical region. We will show in the next subsection, that by making 
use of the integrable structure a± 9 = 0, one can show that the J2 - J3 model is also integrable at tanf? = l/ip. 
We will confirm that the critical point at tan 6* = Ijf (with both couplings positive) is indeed described by the Z3 
parafermions. In addition, we will show that for both couplings negative, the critical theory describing the model is 
that of the Gepner parafermions related to su(3)2 [2S]- This latter integrable point lies in an extended critical region. 
Because there are no relevant operators in the same symmetry sector as the ground state, this whole critical region 
will be described by the same critical theory as the one we found at the integrable point. 

D. Anyonic chains of su{2)k anyons 

In this subsection, we will describe the generalization of the golden chain, where the Fibonacci anyons are replaced 
by more general types of anyons. These more general anyons are of the type which is dubbed su{2)k, where fc is a 
positive integer. For arbitrary k, this anyon theory has fc + 1 types of anyons, which can be labelled in terms of an 
'angular momentum' I, which takes the values I = 0, |, . . . , |. For our present purposes, we are mainly interested in 
the fusion rules of these anyons, and the associated F-symbols, which are necessary to construct the Hamiltonians of 
interacting anyons of this type. 

The fusion rules of the su{2)k anyons are derived from the tensor products of spin representations of SU{2). These 
have to be modified, to take into account that in the anyon model, there is a highest angular momentum. This 
generalization reads as follows. The fusion of two anyons of type ji and is 

min(ji+j2,fc-jl-j2) 

ji X j2 = E , (5) 

i3 = |jl-j2 

where the sum is either over the integers or half-integers. The only difference between the tensor product rules for 
SU{2) spins is the upper bound. In particular, the case of the Fibonacci anyons corresponds to fc = 3. In general, this 
theory has four anyons, / = 0, 1/2, 1,3/2, but because fc = 3 is odd, one can restrict oneself to the integer subset (see 
Ref. [20] for details). This integer subset was written as {1,t} in the previous subsections. 

The model one now considers is the model where the constituent anyons are the I = | anyons of su{2)k. The Hilbert 
space consists of all labelings of the fusion tree in figure 1 but with the t particles replaced by the I = | anyons, and 
at the vertices, the fusion rules in Eq. ([s]) have to be satisfied. The construction of the interaction matrices, both 
for the two-body as well as three-body interactions, is identical to the construction in the Fibonacci case. The only 
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thing which has to be changed is the F-matrix elements. Symbolically, we can write the elements of the resulting 
projection matrices in the same way, 

/p(l) x,_i,i,ixl 
/-pCl/a) x:E.-l,ii,ii+l,a;i+2 _ ^ / , i , i x / i x £i+l / , i x 1/2 , , i , i ^ 



Here, we wrote the projectors onto the spin-1 and spin-1/2 channels for the two- and three-body interactions respec- 
tively. The F-matrices themselves can be obtained from the explicit expressions, which were derived in Ref. [21 . The 
Hamiltonian for general k now reads 

Hj^.j, = Z cos eP^'_ldy,^ + Sin ePa-'body. ' (8) 

i 

where the k dependence is hidden in the detailed form of the projectors, and of course in the Hilbert space itself. 

The phase diagram of the general k model has the same structure as the phase diagram for fc = 3 in Fig. |4] The 
extended critical region around = is described by the minimal model Mk+i,k+2, the fc-critical Ising model. The 
critical phase around = tt is described by the Zk parafermions (we refer to Ref. for a description of this CFT) . 
At angles 9 = 0,tt, this follows from the integrability of the RSOS model |2]. 

We introduce the notation dk = 2cos(^n/{k + 2)) for the quantum dimension of the spin-1/2 anyon of the su{2)k 
anyon model. Below, we will show that the angles given by tan0 = (d^. are special, because we can obtain the 

critical behavior by mapping the model to a new integrable generalization of the RSOS model. For the resulting 9 
in the range < 9 < 7r/2, this integrable point is the transition from the extended critical region to a gapped phase. 
From numerics, it was already obtained that this critical behavior is described by the parafermion theory [3J. 
Below, we show that this indeed follows by exactly solving the model. The opposite point lies within the so-called Z4 
critical region, with low-lying states at the momenta K = 0, tt/2, tt, 37r/2. For this critical region, the numerical results 
were less clear, but our analysis of the integrable 2D classical statistical mechanics model, which we introduce in this 
paper, shows that the conformal field theory description of this phase is in terms of a diagonal coset model, namely 
sM(2)ixsM^2)ixait(2)fc_2 ^ Yqj. k = S, this model reduces to the su(3)2 parafermion CFT. Bordering this extended critical 
phase is the analog of the Majumdar-Ghosh phase around 9 = 3tt/2, which also borders the extended critical phase 
around 9 = 0. The phase transition between the latter two is described by the k + 1-critical Ising model. 



III. CONNECTION WITH 2D STATISTICAL MECHANICS MODELS 



A. Integrability of the Golden chain model 



In the original paper |2] , it was pointed out that the golden chain Hamiltonian can be solved exactly by mapping it 
onto the restricted solid-on-solid model, which was exactly solved by Andrews, Baxter and Forrester [8] by means of the 
corner transfer matrix method In particular, the two-body terms ^'2^body i Hamiltonian at = 0, Hj^^ij^^q = 

-Pi-body i' '^^^ related to generators of the Temperley-Lieb algebra e(i), namely e{i) = dfc(l - -P2^body J- 
The Temperley-Lieb algebra generators satisfy the relations 



e^(i) = dk e{i) 
e{i)e{i ± l)e(i) = e{i) 

[e(z),e(j)] = for \i-j\>2 



(9) 



The action of the Temperley-Lieb generators on the local degrees of freedom labeling the states in the Hilbert space 
can be written, following Pasquier [27] . as 



\^\'')xi-i)j.. "Xi-i,Xi+i\ „ „ 

\ JO,Xi_iJO,x,+ i 



(10) 

(11) 
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Figure 5: The plaquette weight W{li,l2,l3,l4), which in tlie anisotropic limit (it 0) gives rise to the two-body Hamiltonian 



where 5*;^ are the elements of the modular S'-matrix of the su{2)k conformal field theory, which are labeled by 
i,j = 0, 1/2, . . . , k/2, corresponding to the different type of anyons in su{2)k theory. Explicitly, one has 



M k + 2 \ k + 2 ' ^ ' 

We can write down plaquette weights, or _R-matrix elements of the corresponding 2D classical statistical mechanics 
model, in terms of the e(i) and the identity operator 1 as follows 

R^{u)i = iP^^ti + ^l^e(*)r) E [116,, )w{x,.,,x',,x,,,,x,) (13) 
V sm(^) sm(^) > ' ^1 

The subscript i on the i?-matrix labels the plaquette on which it acts, while the argument u is the so-called fugacity, 
and W{xi-i^x'^^Xi+i^Xi) gives the weight of an elementary plaquette shown in Fig. [s] The above i?-matrix can be 
shown to satisfy the Yang-Baxter equation: 



Rj{u)Rj+i{u + v)Rj{v) = Rj+i{v)Rj{u + v)Rj+i{u) , (14) 

by making use of the Temperley-Lieb algebra relations for the e{i) in Eq. ([9|. 
From the i?-matrix, one can construct the 'two-row' transfer matrix (see Fig. [6| 

r(u)=n^2,>i(")n^2,(u) , (15) 

3 j 

which has a role of a discrete time-evolution operator for the corresponding ID quantum system. The time-evolution 
operator acts on a Hilbert space, which is spanned by vectors (..., l2j-i,hj,hj+i,---)- For the Golden chain model, this 
Hilbert space coincides with that given in Fig. [l] We note that in making the connection between the 2D statistical 
mechanics model and the ID quantum Hamiltonian, we have rotated the plaquettes by 45 degrees, or the time runs 
from the south-west to north-east corner of the plaquette. In addition, we consider a 'two-row' transfer matrix, in 
order that in one discrete time step, all heights can evolve. This allows for the possibility to obtain a translationally 
invariant Hamiltonian. The following calculation shows that this is indeed the case. 



hj-i '2j+i 




Figure 6: 'Two-row' transfer matrix, which plays a role of the discrete time-evolution operator for the corresponding ID 
quantum system. The transfer matrix acts on the Hilbert space spanned by vectors {...jhj-ijhj jhj+i.---)- (Imaginary) time 
runs vertically. 
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/ + 1 I 

ai 

I l-l 



l-l I 
ai 

I l + l 







Figure 7: Six height configurations occurring in the RSOS model. 



If Ri{u) satisfies the Yang-Baxter equation, it follows that the T's at different parameters u commute, allowing one 
to construct a Haniiltonian which can be solved exactly. In particular, one writes T(u) = g-"^+o(« ) ^ which gives rise 
to the Hamiltonian 



H ■ 



d\nT{u), 
du 



ln=0 



■E 



1 



dRi{u) I 



Ri{u = 0) du 



lit=0 



Applying this construction to the i?-niatrix defined above, one obtains 



H 



tan(7r/(fc + 2)) Y V2 d 



E 



tan(7r/(fc + 2)) 



y/p(i) _i\ 

Z-i y 2-body,i 2 / 



(16) 



(17) 



which is, up to a •positive scale factor and an overall shift, equal to the golden chain Hamiltonian, -ff ,72=1^73=0 • 

We now focus our attention on the RSOS model we referred to at the beginning of this section. The model consists 
of height variables (simply called 'heights') located at the vertices of the square lattice. The heights can take the 
values / = l,2,...,r-l, where r is an arbitrary integer. We already noted that the heights correspond to the different 
type of anyons, 0, 1/2, . . . , fc/2, where k = r - 2. The connection with the anyon Hamiltonian becomes complete by 
following identification: r = k + 2, and li = 2xi + 1, where li is the value of the height at the vertex i. 

The heights have to satisfy the constraint that they differ by one if they are nearest neighbors. Weights are assigned 
to the different types of plaquettes, which we introduce below. 

This model can be solved for a two-parameter family of weights, namely for the parameter u, and an additional 
parameter p. This parameter p is the parameter which drives a phase transition, located at p = 0. The golden chain 
is related to the RSOS model at this critical point. 

We now briefly describe the weights of the RSOS model in terms of the parameters^ r] (which is related to r, 
see below), u and p. There are six different arrangements of heights around a plaquette, as shown in Fig. [Tj The 
corresponding weights are given in terms of elliptic functions as follows 



ai{u) 



h{2f] - u) 

h{wi + u) 
h{wi) 



h{u) [h{wi.i)h{wu^)fl^ 
h{2r]) h{wi) 
h{wi - u) 
h{wi) 



(18) 



where wi = 2rjl. The function h{u) is given in terms of elliptic theta functions, with argument u and modulus k = m^, 
namely h{u) = H{u)Q{u). The functions H{u) and Q{u) can be expressed in terms of the theta functions 9i and 64, 
in particular, H{u) = (^i{ 2K(m) 'P^ ^^'^ &{u) = ^'4 ( 57^^^ , p) 1 where K{m) is the complete elliptic integral of the first 
kind. The parameter 77 is given in terms of K and r, namely rj = K/r. Furthermore, p can be expressed in terms of m 
as p = exp [-t: K' (m) / K (m)], where K'{m) = K{1 - m). Using the product expansions of the elliptic theta functions, 
one can write 



(19) 



Note that we suppressed the dependence of the weights in Eq. ( 18 1 on p. We refer to chapter 15 of the book ^ , where 



the properties of the elliptic functions used in this paper are analyzed. It was shown in Ref. [8^ that the weights given 
in Eq. ( 18 I satisfy the Yang-Baxter equation for all u and p. 



^ The parameter u is related to v in Ref. [8 by u = ?j - 
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Two phase transitions (in different regimes for u) occur for p = 0, which implies that m = 0. These two critical 
points correspond to the integrable point of the golden chain model with tan6' = {i.e. the angles 9 = 0,tt of the 
J2 - J3 model). Thus, to relate the weights in Eq. (18 1 at the critical point to the plaquette weights obtained from 
the golden chain model, we take the limit limp^o h{ui) / h{u2) = sin(ui)/ sin(u2), where we used that K{0) = 7t/2. 

Using this limit, one finds, at the critical point, that 
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sin(^+») 



Bin(u) [sin((^)sm((^)]V2 



sin( 



(20) 



k+2 



u) 



which are identical to those obtained from Eq. ( 13 1 



( 13 1, when Xi^i + Xi 



The last statement can be verified by considering the explicit form of the i?-matrix in Eq. (|13l 
plaquettes are obtained from the first term in Eq 
terms, with Xi + x'^, which can occur if Xi-i = Xi+i. 

Finally, the plaquettes of type 7 and 5 (which are diagonal, but also have x 
terms in Eq. ( 13 1. 



The a type 



The P plaquettes are the 'off-diagonal' 
Only the second term in Eq. ( 13 1 contributes to these plaquettes. 

- xm) receive contributions from both 



Details of the connection between the critical behavior of the golden chain can be found in the original paper [5]. 
For details of the various phases of the RSOS model, we refer to Refs. and jTU]. Approaching the critical point 
p = from the positive p side, one finds that the observables of the model, such as the height probabilities (see below), 
are given in terms of the tri-critical Ising model for it > 0, and the three-state Potts model for w < 0. This establishes 
the observed critical behavior of the golden chain model [5]. In general, i.e. for arbitrary k, this generalizes to the 
/c-critical Ising model for antiferromagnetic interactions (u > 0), and the Zk parafermions for ferromagnetic (u < 0) 
interactions. 



IV. CONSTRUCTION OF A NEW COMPOSITE HEIGHT MODEL 



A. R-matrix for the composite model 

To identify a new integrable model, which corresponds to different parameter regimes of the generalized golden 
chain model (J2 - J3 model), we will use the ideas put forward in two papers by Ikhlef et al. |17||T8) . In these papers, 
a model closely related to the anyonic chains is studied. Namely, the underlying algebraic structure, the Temperley- 
Lieb algebra, is the same, but a different representation is chosen. We consider the 'anyonic representation' (see Eq. 



(Ill), while Refs. [171 ITS] consider a spin-l/2-type representation. An important consequence of this difference in 
representation is that we were forced to use the corner transfer matrix method solve the model, as opposed to the 
Bethe Ansatz method pS], which was used in Refs. fl7| ITS]. 

Below, we will introduce a new height model, which at its critical point reduces to the anyonic chain at the integrable 
point. However, we will also study the ordered phases of our new height model. 



It was put forward in Ref. |T7], that one can use the i?-matrix in Eq. (13 1 to construct a composite i?-matrix R, 
which also satisfies the Yang-Baxter equation^. The composite i?-matrix one has to consider takes the following form 

Rj{u,(l>) = R2j + l{u - (j))R2j{u)R2j+2{u)R2j + l{u + (f)) , (21) 

where we introduced an additional parameter 0, ranging over < (f) < 7r/2. It can easily be shown that the i?-matrix 



(21 1 satisfies the Yang-Baxter equation. The only ingredient needed to show this is that the original _R-matrix satisfies 



the Yang-Baxter equation itself. In what follows, we focus on the case (j) = 7r/2, but we stress that the i?-matrix in Eq. 



(21 1 satisfies the Yang-Baxter equation for all values of ip. At (j) = 7r/2, the composite i?-matrix leads to particularly 
interesting points of the J2 - J3 Hamiltonian, namely the critical point between the 'Haldane gap' phase and the 
extended AFM critical region. The opposite point (considered in Refs. |17L [TB] in the different representation) lies 
within the Z4 critical region. 



^ This way of constructing a composite i?-matrix reminds of the techniques to construct higher spin, or 'fused' models, see Refs. 1121 1291 
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(22) 



We now describe in some detail how to obtain the Hamiltonian of the J2 - J3 model, by taking the 'anisotropic 
limit'. We start by expanding the composite matrix Ri{u^(f>) explicitly (using the notation 7 = 7r/(/c + 2) and Ci = e(i)), 

sin^{^)Ri{u, <p) = - (sin(7 - u)^ sm{(/) + 7-1*) sin((/) -7 + 1/))! 

- (sin(u) sin(7 - u) sin(0 + 7 - m) sin(0 - 7 + u)){e2i + 624+2) 

+ - sin(u) sin(7 - m)(1 + 2cos(20) - 003(27) - cos(2(7 - u)) - cos(2u))e2,;+i 

+ (sin(u) sin(7 - u) sin(0 - u) sin((/> - 7 + u)){e2ie2i+i + 621+2621+1) 
+ (sin(u) sin(7 - u) sin(0 + u) sin((/) + 7 - M))(e2i+i62i + 62^+162^+2) 

- sin(M)^ sm((j) + 7 - m) sm{(j) - 7 + u)e2i62i+2 
+ sin(M)^ sin(0 - u) sin(0 - 7 + ■ti)e2i62i+2e2i+i 
+ sin(M)^ sin(0 + u) sin(0 + 7 - u)e2i+i62ie2i+2 

- sin(M)^ sin(0 + u) sin(0 - u)62i+ie2j62i+2e2i+i 

We note that the coefficients of the terms 62^62^+162^+2 and 62i+262i+ie2i are zero. 

As we explained in the previous section, one can construct a Hamiltonian related to this i?-matrix via the 'two-row' 
transfer matrix T{u) = Y[j R2j+i{u) Y[j R2jiu) = e""^+°(" \ by taking the anisotropic limit, Eq. ( [T6| . Applying this 
procedure to the composite i?-matrix Ri{u,(j)), one sees that in the 'two-row' transfer matrix, see Fig. [g] one has to 
change R to R. Consequently, it is not obvious that the procedure yields a translationally invariant Hamiltonian. 
Applying the procedure, one obtains 



^^ 2cos(27)-cos(2^)-l ^^ 1 , 
~^ 2sin(7 + 0) sin(7 - 0) sin(7) 



sin(7) sin((/) + 7) sin((/) - 7) 



(sin((/)) cos(7)(ejej+i + e^+ieO + (-l)*cos(^) sin(7)( 



616^+1 6^+16. 



)) 



(23) 



To make the connection with the anyonic chain, we focus on the case (/> = 7r/2, which gives (after dropping the irrelevant 
constant) 



H 



sin(27) 



- C0s(7)(6i + 6j+l) + (6^6^+1 + eiCi+l) 



(24) 



We note that for = 0, we obtain the original golden chain model H = - 6i. In the case that < <j) < 7r/2, the 
term (eie,;+i - e.^+iei) has a non-zero coefficient, which gives rise to a non-hermitian Hamiltonian, and breaks the 
translational invariance. 

To relate the Hamiltonian we just obtained from the composite i?-matrix to that of the J2 - J3 model, we write 
the two- and three-body projectors appearing in the J2 - J3 anyonic Hamiltonian in terms of the Temperely-Lieb 
generators 6^. In general, the projectors of p 'spin-1/2' particles onto the 'spin-p/2' channel, can be written in terms 
of the 6i. Explicitly, for two and three particles, one has (see also [30j [31]) 



2- body,'i 

d(3/2) 

3- body,2 



1 



1 



dl-l 



(25) 



where we remind the reader that dk = 2cos(7r/(fc + 2)). In the anyonic spin chain the three-body interaction was 
written in terms of the projector onto the spin-1/2 channel, which reads 



p(i/2) ^ flfc 

3-body,2 ^2 _ -j^ 



(Ci + 6i+i) 



1 



(CiCi+i + 6i+i6i) 



(26) 



We can now write the J2 - J3 model in two different ways, namely 



H 



J2-J3 



(1/2) 



(27) 
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He = ^cos(6'e)(ei + Ci+i) + sin(6'e)(eiej+i + e^+iCi) 



(28) 



By making use of the projectors in Eq. (25 1, we find the following relation between the angles 6 and 9e 



tan ■ 



{dl - 1) tan 6*6 
2dk{l + dk tan 6*6) 



tan dp 



{dl-l)-2dltanl 
2di. tanO 



(29) 



The angle 6^ for which the model is integrable, 6'e^int, can be read off from Eq. (24), and is given by tan0e,int 
-1/ cos 7 = -21 dk- This corresponds to the angle 6'int in the J2 - J3 model 



tan^i, 



4-1 



(30) 



In particular, in the case u > 0, we find that the corresponding angle 9 is given hy 9 = arctan((d^ - \)ld\) + tt, i.e. 
when both the two- and three-body interactions are ferromagnetic. When m < 0, we have 9 = arctan(((i| - l)/d^), and 
both interactions are antiferromagnetic. 



B. Constructing the composite height model 



As we described in the previous section for the original golden chain model, we have to consider a more general 
two-dimensional height model, in order to obtain the critical behavior of the J2 - J3 model at the integrable points. 
The plaquettes of this new integrable height model, which we introduce below, reduce to the composite i?-matrix 
described above at the critical point p = 0. 

From the construction of the i?-matrix, we know how to construct the plaquette-weights for the composite 2D 
classical statistical mechanics model from those of the RSOS model described in the previous section. The new 
plaquettes consist of four plaquettes in the original RSOS model, as depicted in figure [8] The model again lives on the 
square lattice, where both the vertices as well as the middle of the links have a height variable. Two neighboring heights 
have to differ by one, as in the original model. One can think of these plaquettes as composite plaquettes, whose 
weights depend on the four original plaquettes forming the composite one. The original plaquettes each contribute to 
the weight of the composite plaquette, but two of the plaquettes have an appropriate 'shift', as in the construction of 
the composite i?-matrix (see Eq. (21)). 




E 



+ 



Figure 8: The plaquette weight W(/i, ^2, ^3, '4, ^5, ^e, ^7, ^s) for the composite model in terms of the weights of the RSOS model. 

Using the six different types of plaquettes of the RSOS model, it turns out that one can construct 66 different 
types of composite plaquettes, which are given in Appendix]^ Not all the weights of these plaquettes are different. 
Moreover, they satisfy certain symmetries, in the same way as the original RSOS model. The actual number of different 
plaquettes depends on the parameter r, which determines the number of allowed heights, via Z,j = 1, 2, . . . , r - 1. 

The composite weights W (see Fig. (Isl) in term of the RSOS weights have the following form: 



W^(^2j-lj '2^' '2j+lJ ^2j"+2! ^2j+3i l2j+2ihj+lj2j) 



(31) 



X! ^ l2j+2,l23 + l)W{l2j-l,l'2j,l, l2j)W{l, ^2j+2 i ^2j+3 , ^2j+2) (^2j , ^2j + l ; ^2j+2 : 0> 



where MK*(u) = W{u ± K) are the weights of the shifted plaquettes (see Fig. [s]). We note that each of the plaquettes 
forming the composite plaquette is of type a; , /3; , 7; , (5; , with the appropriate values of I, namely 

W{1 + 1,1,1 - 1,1) = W{1 - 1,1,1 + 1,1) = ai{u) , 
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W{1, l-l,lJ + l) = W{l + l,l,l- 1, 1) = I3i{u) , (32) 
Wil,l + 1,1,1 + l)=^i{u) , 
Wil,l- 1,1,1 -l)=di{u) , 

where the expHcit expressions of these weights, in terms of the parameters u and p, are given in Eq. ([Ts]). 

Taking into account the quasi-periodic properties of h{u), we have chosen the shift (p to equal K. Again, similar to 
the RSOS model, we will be interested in p ^ (critical) limit and note that in this limit, (j) = K = tt/2. Due to the 
symmetry properties of elliptic functions we only need to consider the region 2rj-K<u<2ri + K, which naturally 
breaks into two domains according to the sign of u (cf. Ref. [8]): 



Vi : Q<u<2i^ + K = {2 + r)'q , (33) 
2?2 : 2-q-K={2-r)ri<u<Q . (34) 

The different signs of the fugacity correspond to ferro- (2?i) and antiferromagnetic (2?2) regimes, which in the 
anisotropic limit give rise to the integrable point of the generalized J2 - J3 model, for both signs of the interac- 
tion. 

We have now completely specified our new height model. We employ the corner transfer matrix method, which is 
described in the next section, to solve it. The main interest is to calculate local height probabilities in different domains 
for general p, e.i. away from criticality. It has been observed that these off-critical {p + 0) local height probabilities of 
an integrable lattice model can be mapped to partition functions of the corresponding critical theory {p = 0) in a finite 
box with appropriate boundary conditions [llj . This mapping is realized if one properly relates p, which plays the 
role of temperature, to the finite size L of the critical system. In addition, it has been realized that the local critical 
probabilities can be written using characters of the underling CFT jT2]. Relying on these observations, we identify 
CFTs describing critical theories of the generalized anyon model [p = 0) in subsequent sections. In particular, we 
calculate the off-critical local height probabilities and relate them to characters of a CFT, which governs the critical 
properties of the generalized anyon model (as well as the generalized RSOS model at p = 0). 



V. CORNER TRANSFER MATRIX METHOD AND LOCAL HEIGHT PROBABILITIES 



A. Definition of corner transfer matrices 



To exactly solve the generalized model, we use the corner transfer matrix method in analogy to the solution of 
the RSOS model by Andrews, Baxter and Forrester [H]. Here we give a short account of the method and turn the 
interested reader to literature for more details [HIES- The object of interest is the local height probability Pa, which 
is the probability for a site to have height a. This height probability is given by 

configurations plaquettes 

where product is over all plaquettes (faces) of the lattice and sum runs over all allowed 2D height configurations, and 

m 

{Sa)u' = Kh.a)Y\5{k,l[) . (36) 
1=1 

The size of the system is parametrized by to, which should not be confused with the m related to the modulus m? = k 
of the theta-functions which appeared in section III The meaning of the indices 1, 1' of the matrix S will become clear 
shortly. 

The partition function Z, which is given by 

configurations plaquettes 

can be expressed as 



Z = Tr{ABCD) , 



(38) 
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h = a I2 Ic 



I'o 



d 



Figure 9: Comer transfer matrix A for finite m, constructed using the composite plaquettes. Heights b,c,d,e are the boundary 
heights fixed by a ground state pattern chosen. 



by introducing corner transfer matrices A, B,C, D, corresponding to lower-right, upper-right, upper-left and lower- left 
quadrants of the lattice (see Fig. ([9|, and text below for the precise definition of A,B,C,D). These corner transfer 
matrices are analogs of the row-to-row transfer matrix T, but instead of adding a row to the lattice, they add a whole 
corner. 

Finally, the local height probability can be written as 



Pa = Tr{SaABCD)/Tr{ABCD) . (39) 

We now have to introduce the notion of ground states, which will be used extensively in the subsequent analysis. 
Ground states are those configurations of heights, which maximize the summand (or minimize the "energy") in Eq. 
(37 1. These ground states depend on w in a way that in different domains of u (I?i and 'D2) different ground state 
types exist and they change discontinuously across the boundary between different domains (u = 0). These different 
ground states determine different critical behavior in corresponding domains as we shall show below. 

The corner transfer matrices can be expressed using local plaquette (face) transfer matrices 



{Uj)\,\' - W{l2j-l,l'2j,l'2j + lJ'2j+2^l2j+3il2]+2,hj + l,l2j) Yl HhJi), (40) 

i=l*2j,2j+l,2j+2 
rn 

(^7)1,1' = ^{l'2] + lJ'2j+2J2j+3il2j+2,l2] + l,hj,hj-l,l'2j) HhJi)- (41) 

i=l*2j,2j+l,2j+2 

Here j = 1, (m+ l)/2, with m odd, and Im+i, I'm+ii Un+2, Kn+2i ^-m+s, Kn+3^ ^m+i are boundary heights, which should be 
fixed to proper ground state values corresponding to a considered domain of u. The matrix Uj (Vj) adds a plaquette 
to the lattice in the NE-SW (NW-SE) direction. 
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Using the definitions in Eqs. ( 40|41 ), the corner transfer matrix A can be expressed as 

A = FiF2...F^m^iy2 , (42) 

where 

Fj = U(^rn+l)/2U(m-l)/2---Uj ■ (43) 

The corner transfer matrix A (see Fig. |9| has rows and columns, which are labeled by the values of the boundary 
heights, collected in the vectors 1 = (hjh, ■ ■ ■ Jm) and 1' = {1[,1'2j ■ ■ ■ j^'m)i where l[ = li = a, because these heights 
correspond to the same, central height. In addition, we will use boundary conditions, such that the heights on the 
four outermost diagonals are fixed to be 6, c, d, e, respectively. 

The corner transfer matrices B,C,D are expressed similarly to A by replacing Uj with Vj, Uja,nd Vj" , respectively. 
In general, properties of the corner transfer matrices will depend on symmetries of weights as well as those of ground 
states. In what follows, we will be interested in infinite lattice limit, m ^ oo. 

B. Corner transfer matrices as exponentials. 



From the Yang-Baxter equation for weights (see Eq. (14)) follows a very important property of corner transfer 



matrices. In the limit where the lattice size goes to infinity we can write (symbolically) that 

lim BC= lim T", (44) 

where T is the row-to-row transfer matrix and n is the number of rows, covering a half plane in the limit n -> oo. 

We should note that the above relation is not valid for finite m and n, since different boundary conditions are 
used to calculate left- and right-hand sides of the equation (in fact, even the shapes of the lattices differ). However, 
in the large m, n limit this difference becomes negligible. The Yang-Baxter equation ensures that the row-to-row 
transfer matrices with different fugacities commute and, hence, the product B{u)C{v) depends only on the difference 
u-v (modulus the overall multiplicative factor) . Similar equations that involve other corner transfer matrices can be 
obtained by rotating the lattice in steps of 7r/2. Using these properties one can show that the corner transfer matrices 
have the following form (dropping irrelevant multiplicative factors) [3]: 

A{u) = QiMie""«Q2' , (45) 
B{u) = QsMae^^Qgi , 
C{u) = QsAhe-^^Q-^' , 
D{u) = QiMie'^'^Ql^ , 

where the matrices "H, Qi, Q4, Mi, M4 are independent of u and can be chosen to commute with 5*1, Sr-i- In 
addition, the matrices "H, Afi, M4 are diagonal. 
Using the identity (see Appendix [b| 

A{0) = QiAfiQa' = 1 , (46) 

we immediately see that 

A{u) = Q^e-^^Q-^' , (47) 

which implies that the diagonal form of the corner transfer matrix A can be written as an exponential. From Eqs. 
( 39|45 1 we immediately see that 



Pa = Tr{SaMiM2MM/Tr{AhM2M3A'U) . (48) 



To calculate M1M2M3M4 we need to consider different domains of u separately to find different limiting properties 
of the corner transfer matrices. Here we summarize the properties of the corner matrices near the boundaries of 
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domains Vi and 2?2 (see Appendix [B] for details), which are required for the calculation. In the domain Vi, in the 
M -> limit we find (upon dropping an irrelevant multiplicative factor) 

A{u = 0) = C{u = 0) = t, (49) 

while in the u ^ (2 + r)r; limit, we have 

B{u={2 + r)r]) = D{u={2 + r)rj) = Vi, (50) 

where we have defined 

{V,),,v = [h{2r^h)Y''S{l,V). (51) 
Using the above limits for corner transfer matrices we can write that 

A{u = 0)B{u = (2 + r)r])C(u = 0)D{u = (2 + r)7j) = V\. (52) 



Substituting Eq. (45 1 in the above equation we get the following result for the product of the matrices M1M2M3M 



4, 



which appear in the expression for the height probability in Eq. ( 48 1 

In the domain 2)2, we can show that 

A{u = 0)B{u = (2 - r)Tj)C{u = 0)D{u = (2 - r)?]) = V\ (54) 

and 

MiM^M^iMiC^^^-''"''^'^ = V\. (55) 



Taking into account relations in Eqs. (53) and ( |55[ ), the local height probability can be written as 

Pa = Tr{SaVle-^''^^)/Tr{Vle-^*'''^), (56) 

where 

\2-r ueV2 



C. Diagonal form of the corner transfer matrix A. 

What is left to do is to find the diagonal form of the corner transfer matrix A and the matrix H, which is the most 



involved part of the calculation. The corner transfer matrices in Eq. (45) should satisfy quasi-periodic conditions 
as do the weights, with the period 2iK' , which implies that elements of H are integer multiples of ir/K' , H\,v = 
7rN{\)5{l,l')/K' . Similar to the solution of the RSOS models fS], we assume that Ti does not change discontinuously 
with p, which implies that the integer function A^(l) is independent of p and we can derive it in a limit where the 
composite weights assume a simple form. We can show that (see Appendix [C| in the p ^ 1 limit the weights of the 
composite model take a particularly simple form and the corner transfer matrix A can be readily diagonalized. After 
a fairly lengthy calculation, involving the 'conjugate-modulus' (or modular) transformation, the dust settles, and one 
finds the diagonal form of A 

Au' = [e-"«],,r = gr>^^'^^''5(U'), (58) 

where 

«; = e-2-"/^' = ^i;(2;i-r)^/(i6.) (59) 

and 

{m+l)/2 |j _ j I 

iV(l)/2E0(l)= ^ ^j l2,.3- '^-'K s,,^_,,,^^,Si,^^,,,^^A,^,,^^^, (60) 
j=i 

where sum over j is performed along a line in the 2D lattice (see Fig. |9|. Each term in this sum corresponds to 
the weight of the j*^ plaquette (counted from the central site), times j, which is the number of plaquettes on the 
J*'' diagonal. In the limit p = 1, these plaquettes all have the same form, because A is diagonal then, as shown in 
Appendix O 
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D. Local height probabihty. 

We can now collect the results, and obtain the local height probabilities Pa- Substituting the diagonal form of A, 



Eq. (|58| in Eq. (|56f , we find that 

Pa = S'^VaXm{a;b,c,d,e]x*). (61) 

Here we adopted the following definitions: 

XUa;b,c,d,e;q)= ^ (62) 

^^ = ^(2-0(2a-O^/(16r)^(^a^^.^^ (63) 

S = Y,VaXrn{a;b,c,d,e;x*), (64) 

a 

X = e^^'^"/-^' . (65) 



The function (f){V) was given in Eq. (60 1 and the function E{z,x) appears in Jacobi's triple product identity, 



E{z,x) = n(l-^'"'2)(l-x"z-i)(l-x") . (66) 

ri=l 

Furthermore, li = a and heights b, c, d, e (see Fig. [9| are the boundary heights, which are to be fixed to the values of 
a ground state in the domain of u under consideration. We see from Eq. (62 1 that the partition function as well as 



the local height probability is expressed as a sum over ID height configurations, in contrast to 2D configurations in 



the original formulation (see Eq. (37l). This property is the consequence of integrability of the model and greatly 
simplifies calculations. 

VI. PHASES OF THE COMPOSITE HEIGHT MODEL 

In this section, we will make a start with the exploration of the phase diagram of the composite height model, 
which will be focussed on the regions which are related to the anyonic quantum chain. We will give a more detailed 
description of the various phases of the model in a forthcoming publication. 

The phase diagram of the composite height model bears resemblance to the phase diagram of the original RSOS 
model. We will consider the phase diagram as a function of the parameters u and p, where the parameter u is related 
to the anisotropy of the model. Only the sign of this parameter will be relevant. The parameter p drives a phase 
transition between different ordered and disordered phases, as was the case in the RSOS model. We consider the 
regime < p < 1, which is the one relevant for our purposes. 

We start by analyzing the ordered phases by taking either u > 0, which we call regime III, or w < 0, called regime 
II. The naming of the regimes follows the nomenclature of ABF. We use the results of the previous section, where we 
calculated the height probabilities in terms of the functions 

X„,{a;b,c,d,e;q)= ^ /^'^ , (67) 

where the vector 1 has m + 4 components, 1 = (a, Z2, • ■ • , Im, b, c, d, e), which satisfy the constraint li = li^i ± 1, implying 
that a + e = mod 2. We repeat the function 0(1) for convenience, 

(m+l)/2 17 _ J I 
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1232 123 21 

Figure 10: Ground state pattern for u> 0, which is characterized by the pattern in the dashed box. 
For < p < 1; the height probabilities are proportional to X,„(a; b, c, d, e; x*), with i = 2 + r for m > 0, and t = 2 - r for 

M<0. 

To find the ground states, we analyze the function 0(1). The ground states are those configurations which contribute 
maximally to the partition function. In domain 2?i, which has u > 0, and t = 2 + r > 0, the ground states are given by 
those configurations which minimize the function 0(1). For domain 2?2, with m < and t = 2 - r < 0, the function 0(1) 
should be maximized instead. As long as < p < 1, one finds that the arguments about the ground states go through, 
because x < 1. At the critical point (p 0) we have a; ^ 1, hence the argument fails and all height configurations 
contribute equally. 

The model is critical when p ^ 0, and we will study the full height probabilities Pa, which give the probability that 
the central height takes the value a, depending on the boundary heights {b, c, d, e) (see Fig. [9]). We will evaluate these 
height probabilities, in the case that the boundary heights are such that they belong to a ground state pattern. 

A. Ground states for u> (domain Di) 



Let us start by analyzing the case u > 0, and minimize the function 0(1). In this case, as many plaquettes as possible 
should give zero contribution to 0(1). This can be achieved in the following way. First of all, one should have that 
li = such that the first term within the parenthesis in Eq. (68 1 is zero. There are now two different ways in which 
one can avoid a contribution from the second term. First, one can set = ± 2 and li+i = = ± 1, where the 
sign in both equations should be the same. The other possibility is = h, while Z^+i = li ± \ and ^^+3 = liT 1, where 
the signs in the last two equations have to be opposite. 

In particular, the vector 1 describing the ground states forO<p<l,w>0 takes the form (we show the case r = 7) 

(1,2,3,2,1,2,3,2,1,...) (2,3,4,3,2,3,4,3,2,...) (3,4,5,4,3,4,5,4,3,...) (4,5,6,5,4,5,6,5,4,...) 
(3,2,1,2,3,2,1,2,3,...) (4,3,2,3,4,3,2,3,4,...) (5,4,3,4,5,4,3,4,5,...) (6,5,4,5,6,5,4,5,6,...) 

for the ground states of the first type. The patterns of the second type are mere translations of the patterns of the 
first type, and are given by the vectors 1 of the form 

(2, 1, 2, 3, 2, 1, 2, 3, 2, . . .) (3, 2, 3, 4, 3, 2, 3, 4, 3, . . .) (4, 3, 4, 5, 4, 3, 4, 5, 4, . . .) (5, 4, 5, 6, 5, 4, 5, 6, 5, . . .) 
(2, 3,2,1, 2, 3, 2, 1,2,.. .) (3, 4, 3, 2, 3, 4, 3, 2, 3, . . .) (4, 5, 4, 3, 4, 5, 4, 3, 4, . . .) (5, 6, 5, 4, 5, 6, 5, 4, 5, . . .) 

In the limit p = 1, the corner transfer matrix A is diagonal, which gives rise to ground states patterns which are invariant 
under translation along the North-East to South- West (NE-SW) diagonal. A particular ground state pattern for u > 
is displayed in figure 10 All the ground states are of the form (^ - 1, Z, Z + 1, /, / - 1, /, / + 1, . . .), or translations of this 
pattern. 
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Figure 11: Ground state pattern for u< 0, which is formed characterized by the pattern in the dashed box. 

To count the number of different ground states, we note that the ground state patterns are specified by three 
consecutive integers. Because the heights can take the values 1,2, ...,r - 1, there are r - 3 possible consecutive 
integers. By translation, each of these sets of consecutive integers gives rise to four different ground states, for a total 
of 4(r - 3) ground states. 

We will now count the number of different height probabilities, Pa{b, c, d, e), where the boundary condition (&, c, d, e) 
corresponds to a ground state pattern. First, we note that a, the height of the central site, and e have to have the 
same parity, a + e = mod 2. For r odd, there are (r - l)/2 odd valued heights, and (r - l)/2 even valued heights. 
In both cases, a can take (r - l)/2 values, giving 2(r - l)(r - 3) different height probabilities. In the case that r is 
even, there are r/2 odd valued heights, and r/2 - 1 even valued heights. Out of the 4(r - 3) ground state patterns, 
2(r - 3) have e even, and 2(r - 3) have e odd. Hence, also for r even, the number of height probabilities to consider 
is 2(r- l)(r-3). These height probabilities are given in terms of the functions Xm{cL]b,c,d,e;q), which satisfies the 
relation (which follows from the symmetry properties of the plaquette weights) 



Xm{a;b,c,d, e;q) = Xm{r - a;r - b,r - c,r - d,r - e;q) 



(69) 



This halves the number of independent height probabilities, which is thus given by (r-l)(r-3), or in terms of 
k = r-2,hy -1. 



B. Ground states for u <0 (domain 1)2) 



We now consider the ground states for u < 0, and look for configurations which maximize the function (f>{l). 
Because li and maximally differ by four, both terms within the parentheses in Eq. (68 1 can maximally contribute 



1. However, in order that the first term contributes for every plaquette, the heights k would have to steadily increase 
or decrease, which is impossible, because the values the k can take lie in the range li = 1, 2, . . . , r - 1. The second term 
inside of the parenthesis in Eq. (68 1 can be 1 for all plaquettes. The only requirement is that hj-i = hj+i = hj+3 and 



t2j+2 



hj-i ± 1- For p = 1, the ground states again are invariant under translations over the NE-SW diagonal. 



thus, for M < 0, the ground states are given by configurations for which the heights stay as constant as possible, i.e. 
they take the form (/, I + 1,1,1 + 1,1, . . .), as depicted in Fig. 11 for a typical example. 

The number of ground state patterns is given by the number of consecutive pairs (both increasing and decreasing) , 
i.e. 2(r - 2). We again need the number of height probabilities we have to consider, or the number of functions 
X,n{a;b,c,d,e;q). We only have to specify (a;d,e), because d and e fix the values of b and c in the ground states. 
There are 2(r-2) pairs {d,e). It turns out that the number of height probabilities is given by (r-l)(r-2), irrespective 
of whether r is even or odd. Thus, there are (r-l)(r-2)/2 = k{k + l)/2 independent height probabilities, because of 
the relation in Eq. ( 69 1 . 
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C. Criticality at p = 



We found the ground states for u > and u < deep in the ordered regime, namely for < p < 1, by minimizing 
or maximizing the function 4>(l), which appears in the height probabiHties, as well as the partition function. At 
the critical point, for p = as in the original RSOS model, all configurations contribute, and one has to do a more 
careful study of the model. In the remainder of this paper, we focus on the full form of the height probabilities, which 
have a close connection with conformal field theory characters, which allows us to unambiguously identify the critical 
behavior of the model, for both regimes m > and u < 0. Because we already established the connection between the 
composite height model and the anyonic chains, we thereby also explain the critical behavior of the anyonic chains at 
the integrable point. A more detailed study of the critical behavior of the composite height model will be left for a 
future publication. 



VII. EVALUATION OF THE HEIGHT PROBABILITIES 



We now turn our attention to the evaluation of the height probabilities. We have shown that the probability that 
the central height takes the value a depends on the boundary condition, which we specified by fixing the boundary 
heights {b,c,d,e). The height probabilities are governed by the functions Am(a; 6, c, c?, e; q) given in Eqs. (67 1 and 
(68 1. We are ultimately interested in the behavior of these functions in the limit m ->■ oo. In that limit, we can let 
the boundary heights correspond to the ground state patterns, which extremize the function (I>{1). 

We therefore consider the functions Xm{o,; b, c, d, e; g), for all possible values of (a; 5, c, d, e) such that the boundary 
heights {b,c,d,e) are part of a ground state pattern. In the limit m ^ oo, these functions will receive contributions 
from all possible configurations specified by the vectors 1. For finite m, these functions are finite polynomials in 
q (or q^/^ times such polynomials). The minimal power of q which can occur is zero, while the maximal power is 
(m + l)(m + 3)/8. We note that these extremal values are not obtained for all choices of (a; b, c, d, e). 

In the following subsections, we will provide explicit expressions for the functions X„i{a]b,c,d,e;q), in the case 
r = 5. We did not yet obtain explicit expressions for r > 5, but by analyzing the functions, we unambiguously 
identified them as the characters of certain conformal field theories. In particular, we checked extensively that in the 
limit of m ^ oo, the functions Xm{q) tend to affine Lie algebra branching functions, or in other words, characters 
of the various primary fields in certain coset models. This is precisely the connection we are after, because these 
coset models describe the behavior of the model at the critical point, and hence the critical behavior of the anyonic 
quantum chains. Before delving into the details, we will first state the results here. 

For u > 0, we find that the functions Xm(a; 6, c, cZ, e; tend to the characters of the coset ^ 
where k = r - 2. In the case that k = 3, this coset is equivalent to a Gepner parafermion theory j5S] based on 
su(3)2, explicitly y^(^i^^^{i)i2 ' '^^'^ characters of this theory were considered in Ref. |33j in the context of a particular 
non-abelian spin-singlet quantum Hall state [34 . The finitizations of the characters considered in Ref. [33j precisely 
correspond to the functions Xrn,{q) we obtain from the integrable model we introduced in this paper. 

For M < 0, one needs to consider the part of the functions X^iq) with the highest powers of q, as explained below. 
In particular, the functions g("+i)(™+3)/8^^(g-i) 

are the functions one needs to consider in the limit m ^ oo. We 
found that they precisely correspond to the characters of the Zk parafermion theory [26 (which is for instance given 
by the coset ^^y^)- The characters of this theory can be found in Ref. [35 . 

For r = 5, we have the interesting situation that the finitization of the characters of the su{3)2 Gepner parafermions, 
i.e. the functions Xm{q), also can be considered as finitized characters of the parafermions. Because both 
theories contain Fibonacci particles, one could say that 'Fibonacci meets Fibonacci', in the same spirit as 'Ising meets 
Fibonacci' [5], which establishes a connection between the theories su(2)2 and su(3)2, which have Ising and Fibonacci 
type fusion rules. 

In the following two subsections, we provide explicit expressions for the functions Xm{a,;b, c, d, e] q) , for finite m in 
the case, r = 5. These expressions are such that the limit m oo can be taken explicitly. 

The identification of the functions A"m(a; 6, c, d, e; q) for all values of r and (a; b, c, d, e) in terms of the CFT characters 
is the subject of Appendix |D] which deals with both cases u > and u < 0. 



A. Probability amplitudes for u> 

In this subsection, we give the explicit form of the functions A„i(a; b, c, d, e, q), for those values of (a; 6, c, d, e) which 
correspond to the ground states for u> (for r = 5). We did not yet prove these results, but we expect that a proof 
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along the lines of the original paper [S] is feasible. Such a proof involves the recursion relations for the functions 
Xm{a', b, c, d, e, g), which we give in Appendix |D] Eq. ( |dT| ). 

We start by giving the expressions for finite m, and take the limit m ^ oo afterwards. We introduce the function 
('?)m = OfcliCl ~ for positive integers m, and (q)o = 1, which appears in the definition of the q-binomials, or 
Gaussian polynomials 



m 


•1 


n 





(9).(<Z)™-. 





for < n < m 
otherwise 



(70) 



The precise form of the following function we introduce was inspired by the (finitized) character of the Gepner 
parafermions associated with su(3)2 (see Ref. [33J), as well as the characters for the original RSOS model (see, for 
instance. Ref. [13 ). In particular, we introduce 



y{m;l2,l3,l4;q) = Y.' 



a.b>0 



2 

a 



2 

b 



(71) 



where we assume that m is an integer, hjs = 1,2 and I4 = 1, 2, 3, 4. The prime on the sum indicates the constraints that 
the argument of the g-binomials have to be (non-negative) integers. We labelled the function with hjsji, because 
of the connection with conformal field theory characters, which we will describe below (the label li can always be 
chosen as li = 1). 

We then have the following results 



X™(l;2,l,2,3;g 
X™(3;2,l,2,3;g 
X™(2;3,2,3,4;g 
X„(4;3,2,3,4;g 
X™(2;l,2,3,2;g 
X™(4;l,2,3,2;g 
X„.(l;2,3,4,3;g 
X™(3;2,3,4,3;g 
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yi 
yi 
yi 
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(72) 



, TO - 3 



1+1 , TO - 3 



2 

I - 

We can now rather easily take the limit to ^ 00, by using the result lim„^oo 



1,1, i;?) 
;l,l,3;g) 



We will assume that to is 



of the form 4p + 3, with p integer (the case m = Ap+1 is only slightly different). By taking this limit, we find that 
limX4p+3(l;2,l,2,3;<7) = 4"(''')^(g) 



lim X4p+3(4;3,2,3,4;<7) = lim X4p+3(4; 1, 2, 3, 2; g) = lim ^4^+3(1; 2, 3, 4, 3; g) = c'^^^^^q) 

p—yoo p—yoo p—^-oo ^ 

lim X4p+3(2;3,2,3,4;g) = lim X4p+3(2; 1, 2, 3, 2; g) = lim X4p+3(3; 2, 3, 4, 3; g) = g-wc™(3)2(^) 

p^oo p^oo p—>-oo 

(3)2/ 



(73) 



lim X4p+3(3; 2, 1, 2, 3; g) = q-'od;''^^)^ (g) 

Here, the functions c|"^'^''^(g) denote the characters of the su(3)2 parafermion theory. This theory has eight fields, 
including the identity field 1, with the character c™*''^^^(g), and three parafermions, ipi, ip2 and 'ipi2, which have 
identical characters c^"^^^^(g). The remaining four fields are three 'spin' fields cti, (T2 and ai2, which have identical 

characters Q^^^^^{q). Finally there is the field p, whose character we denote by Cp"*''^''^ (g). 

In the above, the characters of the fields are 'normalized' such that the first term in the expansion is q^, where h 
is the conformal dimension of the field under consideration. These conformal dimensions are given by /ii = 0, /i^ = |, 
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h„ = j^, and hp = |. Combining all the factors, the expressions for the limit of Xm{q) are series expansions with 
integer or half integer powers of q. 

For our present proposes, it is best to view this CFT as the following coset model: ^ ^pj-^g fields in 

this coset model are labelled by the labels of the constituent factors, $^i''2''3^ where lijhjs correspond to the factors 
su(2)i, and I4 corresponds to the factor 514(2)3. These labels have to satisfy the constraint + Z2 + ^3 + '4 = mod 2, 
which we use to set li = 0, and consider the fields with Z2 + '3 + ^4 = mod 2. 
For completeness, we explicitly give the labels of the parafermion fields 

^0,0.0 -, ^0,1,1 ; ^0,1-0 ; ^0,0,1 ; 

^0,0.0 ffvO.l.l ^0.1,0 ,^0,0,1 

For r > 5, we did not yet obtain closed expressions for the functions Xm{a;b,c,d,e;q). However, one can obtain 
expansions to high order, by making use of the recursion relations satisfied by the Xm,{q), Eq. (Dll. These high order 
expansions can be compared to the branching functions (or characters) of various coset model. In doing so, we have 
established that the Xm{a;b,c,d,e;q;r) are (in the limit m 00) the characters of the cosets 

(we remind that k = r-2). In Appendix [d| we will give the relation between the values (a;b,c,d,e) and the labels of 
the coset fields. 



B. Probability amplitudes for u <0 

The functions Xm{q), with the boundary heights {b,c,d,e) given by the ground state patterns for m < 0, can also 
be expressed in terms of the function y(m; Z2, ^3j '4) which was introduced in the previous subsection. We first state 
these results, and subsequently take the limit m ^ 00, in order to identify the critical theory describing the critical 
behavior of the anyonic quantum spin chain. In particular, for r = 5 

X„(l; 2, 1, 2, 1; g) = (7"^"y(^, 1, 1, 1) 
X„,(l;2,3,2,3;<7) = <z'^2/(^^,1,2,1) 
X„,(l;4,3, 4, 3;g) = (z"^ 2/(^^,2,2,1) 
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(74) 

X„(2; 1, 2, 1, 2; <?) = q^2/(^^, 2, 2, 2) 
X„,(2; 3, 2, 3, 2; g) = (7"^2/( 1, 2, 3) 

X„(2; 3, 4, 3, 4; g) = q"^2/(^^, 1, 1, 3) 

The ground states correspond to the highest possible powers of q, so to make the identification with the conformal field 
theory, we we will have to make the substitution q ^ q ■, and multiply by g s , which is the maximal power of 
q which occurs in the functions Xm{q)- This will make sure that the function corresponding to the vacuum character 
is a polynomial in q, starting with 1. We find that the functions thus obtained tend to the characters of the su(2)3 
parafermions in the limit to 00. The sm(2)3 parafermion theory contains six fields, the identity, with character 
^su(2)3^ two parafermion fields ipi and ■02, with the character c™^^^^, two spin fields ai and (T2, with characters c^"^^^^. 



and finally the field e, with the character q"^^''^. The scaling dimensions are hi = 0, = ^, ha- = and /i^ = |. We 



4, 

finally the field e, with the character q"^^''^. The scaling 

find 

lim g^'^^^X2,,i(l;2,l,2,l;g-i) =cf ')-^(g) (75) 

p—yoo 

lim q 3 X2p+i(l;2,3,2,3;g-i) = lim q ^ X2p+i(l; 4, 3, 4, 3; g'^) = q-«c7'''((?) (76) 

(p+l)(p + 2) _^ „ . ^ _1^ (p+l)(p + 2) 



lim q'-^^X^M^; 1, 2, 1, 2; q-') = lim q^^^ X^M^; 3, 2, 3, 2; q-') = q'^^C^'^^iq) (77) 

(p+l)(p+2) T 1 

lim g^^X2p,i(2;3,4,3,4;g-i) =<zTocr<'^^K'?) (78) 



The identification of the functions X,n{q ^) for w < in the limit to -* 00 is given in Appendix [d| in the case k > 3 
(r > 5). They correspond to the characters of the Zk parafermions |35| . 
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VIII. CONCLUSIONS AND OUTLOOK 



We have introduced a two-dimensional classical statistical mechanics model the critical properties of which corre- 
spond to the integrable points of a chain of su{2)k anyons with competing two- and tree-body interactions ( J2 - J3 
model, given in Eq. p4|). The 2D classical model is a composite height model which is a generalization of the 
restricted solid-on-solid model solved by Andrews, Baxter and Forrester |5] by means of the corner transfer matrix 
method. We have also used the CTM method and have found that, similar to the RSOS model, there are four different 
regimes with two critical points at p = 0. We have studied two new integrable critical points of the anyonic chain at 
tan6'i„t = - which correspond to the p +0, u ^ ±0 limits of the composite height model. 

For ferromagnetic interactions (u > 0, 9 = arctan(((i^ - + tt, J2, J3 < 0) the critical point is described by 

the ^'"(^)i''^"(^)i''^"(^)'=-2 coset conformal field theory. This critical behavior actually describes an extended critical 

region around the integrable point. For antiferromagnetic interactions (u < 0, 9 = arctan(((i^ - l)/d1), J2, J3 > 0) the 
behavior is that of the Zk parafermions. This critical point constitutes the boundary between a gapped phase, and 
an extended critical region, which is described, in general, by the Ai{k + 1, fc + 2) minimal model [2]. 

These CFT identifications stem from the observations that functions X^iq), which define the local height prob- 
abilities are given by characters of the corresponding conformal field theory. The integrable properties of quantum 
ID and classical 2D models are fully defined by the Temperley-Lieb algebra relations in Eq. (|9|, which suggest that 
the Hamiltonian in Eq. (24 1 can be exactly solved in different representations of TL algebra. A particular phys- 
ical interpretation of the model depends on the particular representation chosen to solve it. In our case, we have 
used the representation of su{2)k anyons, which straightforwardly maps onto the composite RSOS model. In this 
representation, an interpretation as a chain of interacting anyons is straightforward. 

A different interpretation as a Q-state Potts model (with s/Q = 2cos( ^j^) = dk, where in this case, k is considered to 
be a continuous parameter) or a six- vertex model was put forward by Ikhlef et. al. |17l I18| , which naturally admits the 
so called 'loop representation'. The J2 - J3 model at the integrable point, 9 = arctan((d^ - + tt, has been exactly 

solved in spin-1/2 representation of the Uq{SU{2)) quantum algebra |17II18| . In this representation, in contrast to the 
'anyon representation', the Hilbert space has a tensor product structure and the J2 - J3 model admits the solution by 
the algebraic Bethe Ansatz method. We should note that in spin or loop representations the ID quantum Hamiltonian 
is non-Hermitian and lacks obvious physical interpretation. Despite this non-Hermiticity, it has been conjectured that 
the part of the spectrum which scales as 1/L is real. In contrast, in the 'anyon representation', which is adopted in 
this paper, the Hamiltonian is Hermitian (real symmetric) and has a physical interpretation as a local interaction 
between anyons. Despite the differences, the critical properties a,t 9 = arctan((d| - l)/d^) + tt found by us are similar 
to those found by Ikhlef et. al. |17l I18j. Namely, they have the same central charge, at least for u positive. The 
precise connection between these models is interesting, and requires more detailed research. We note that a similar 
loop model, related to the anyon model with u < {9 = arctan(((i^ - l)/d^) ), has been studied in Ref. [36 . It was 
conjectured and several arguments were given that the critical behavior of this loop model (in an appropriate limit) 
is described by Zk parafermions. 

The anyonic chains have been generalized in different ways. First, to different type of unitary anyons [4 . It would 
be very interesting to investigate to what extend the methods of 'fused models', as studied in Refs. [T^ and [22] can 
be applied to the anyonic chains of Ref. Chains of non-unitary anyons were also considered [37]. This leads in 
general to non-Hermitian Hamiltonians, which nevertheless have a real spectrum. It should be possible to generalize 
the models considered in Ref. [38J to composite versions, and make a connection with a generalization of the chains 
considered in Ref. [37 . 

In this paper, we have described only half of the phase diagram of the composite RSOS model, positive p regime 
(0 < p < 1). The corner transfer matrix method allows to study the negative p regime (-1 < p < 0) also, which will be 
addressed in a subsequent publication. We have shown that the positive p regime has the interpretation in terms of 
ID interacting anyons. The ID quantum mechanical interpretation of the negative p regime is yet unknown and also 
requires further investigation. 

Last, but not least, we would like to point out that the study of the RSOS model by Andrews, Baxter and Forrester 
gave rise to an interesting set of Rogers-Ramanujan-type identities. Further study of the composite height model in 
light of these identities will be most interesting. 
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Appendix A: The plaquettes of the composite-RSOS model 



The plaquettes of the composite model are obtained by stacking four plaquettes of the original RSOS model (see 
Fig. [7]), and giving them the appropriate weights. In the process of stacking the plaquettes, one has to sum over the 
internal height, which in some cases can take two possible values (see Fig. [s]). 

In the following, we will display the possible composite plaquette types, and give the associated weights in terms of 
the weights of the RSOS model. In this case, the composite plaquette weights have the same symmetry as the original 
model (with (/>= 0), namely the weights are unchanged under exchange of the North- West (NW) and South-East (SE) 
corners, as well as under exchange of the North-East (NE) and South- West (SW) corners. This last property is lost, 
if + 0, X. 

We start by giving the composite plaquettes for which the internal height is fixed by the boundary heights. In case 
the plaquette is not symmetric, we will indicate the amount of plaquettes which can be obtained from the given by 
taking the appropriate mirror image. These mirrored plaquettes have the same weight as the displayed plaquette. 

The weights of the composite plaquettes are given in terms of the weights of the original RSOS model, given in 



Eq. (18). Let CK'") be any weight of the RSOS model, i.e. C = ct^P^li^- Then, Q = Ci{u), = 0('" + ^) 
Cr ~ CzT^ ~ K). With this notation, we can give the weights of the composite plaquettes for ip= K: 
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We now focus on the composite plaquettes for which the internal height is not fixed by those of the boundary in 
general, and hence takes two different values. The composite weight is a sum of two terms, the first term correspond 
to the internal height being I + 1, the second one to the internal height / - 1: 
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We find that the total number of possible composite plaquette types is 66. The total number of plaquettes depends 
on the value of r = A; + 2, as is the case in the RSOS model. 

Appendix B: Properties of corner transfer matrices in different domains 

To find out the properties of corner transfer matrices near the boundaries of different domains we need to calculate 
the weights of the composite model in those limits. 

1. w -> limit 

Using the properties of the elliptic theta functions it is straightforward to show that the weights of the RSOS model 
in the u ^ limit are 



q;j(m = 0) = 1, 

A(w = o) = 0, 

7Kw = 0) = 1, 
5i{u = Q) = 1. 

The RSOS wights with shifted fugacities, (j) = K = r^r, have the form 
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St{u = 0) = -S[{u = 0)=j[{u = 0) 



h{2'ql - rjr) 
h{2r]l) 



We are able to show that in this Hmit all ofF-diagonal (NE-SW asymmetric) weights of the composite model are zero 
and only diagonal (NE-SW symmetric) weights survive. This implies that the corner transfer matrices A and C are 
diagonal in this limit. Of the 66 different type of plaquettes, 50 are off-diagonal. Of these 50 type of plaquettes, 42 
are trivially zero in the limit u ->■ 0, because they contain a factor I3i{u = 0) = 0. The remaining plaquettes are those 
in Eqs. (All I and (A14l, and they are zero because of the relations between the weights given in Eq. (B2|. As an 



lave for the plaquettes on left hand side of Eq. (Alll: P^^^ Sf^^ + S^ fii af = Pi^i +Sil3i = P^Si -5i = 0. 



example, we 

The weights of all the diagonal plaquettes turn out to be the same. For the weights which only have one contributing 
term, we immediately find Wi = • The weights of the diagonal plaquettes which consist of two terms read 

W2 = -( '^^^2r)^^^ ) ~ ( h(27j)li(2r)i) ) h{27]l - 2ri)h{2ril + 2r]) . By making use of the following identity for elliptic functions 
(see, for instance. Chapter 15 of Ref . [9 ) 

h^{2rj - r]r)h^{2r]l) = h^{2ri)h^{2r]l - rjr) + {'qr)h{2i^{l - l))h{2r]{l + 1)) , 

we find that W2 = Wi. Thus, the general weight of the composite model reads, in the limit u ^ 



fh{2ri-jr)Y 
\ h{2rj) I 



Using the definition of Uj , Eq. ( 40 1 , we can show that 



such that Uj is a diagonal matrix. From the above and the Eq. (42 1 it follows that 

A{u = 0) = C{u = 0) = t, 
where we have dropped the irrelevant multiplicative factor, which only depends on r. 



(B3) 



(B4) 



(B5) 



(B6) 



2. M ^ (2 + r)r] limit 

We show that in the limit u ^ (2 + r)ri, the corner transfer matrices B and D are diagonal. For unshifted weights 
of the RSOS model in the u ^ {2 + r)ri limit we get 
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The shifted weights take the form 
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Because a^{u = (2 + r)ri) = aj{u = {2 + r)r]) = 0, most of the NW-SE asymmetric weights are zero. The remaining 
NW-SE asymmetric weights can be shown to be zero, by making use of the properties Eqs. (B7l, (B8l. This shows 
that for u ^ (2 + r)ri, the corner transfer matrices B and D are diagonal. 

As was the case for u = 0, the plaquettes which contribute in the Hmit u ^ (2 + r)r] fall in two classes, the ones with 
one term and those with two terms. Again, the plaquette weights of these two classes can be shown to give rise to 
the same weights, by making use of the elhptic function relation in Eq. (B3l. The final form of the weights is slightly 



more complicated than in the case u = 0, and does in fact depend on to the heights 
W{h,l2, hM^ h, hM){u = (2 + r)ri) 



[h{27^h)h{2T^h)fl^ |' fe(277-7?r) j ^ 



h{2T^h) 



\ h{2rj) 



Using the definition of Vj , Eq. ( 41 1 , we can show that 



Vj{u= {2 + r)rj) 



where we introduced 



{V,),,y = [h{2^l,)Y"S{l,l'). 



From the above and the definition of the corner transfer matrices, Eq. (42 1, it follows that 



B{u = (2 + r)rj) = D{u = (2 + r)??) = Vi, 
where we again dropped the irrelevant (^i independent) multiplicative factor. 
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3. M ^ (2 - r)r] limit 

In this limit all the weights of the RSOS model just change sign compared to the u ^ {2 + r)rj limit. Hence, the 
weights of the composite model are unchanged and we find 



B{u = (2 - r)??) = D{u = (2 - r)??) = Vi. 
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Appendix C: Weights of the composite model in the p ^ 1 limit 

To derive the integer function we consider the limit p 1 where weights of the composite model, and hence 

the corner transfer matrix A, become diagonal. For < p < 1, employing the conjugate modulus transformation, the 
function h{u) can be written in the following way [8 : 



h{u) = rexp 



T:{u-Ky 
2KK' 



E{e 
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and the function E{z,x) is Jacobi's triple product 
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Then the weights of the RSOS model have the form 
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where 
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and v is Si constant independent of I. 

In what follows, we show that in the limit of p ^ 1 and w ^ 1 the weights of the composite model become diagonal 
(similar to the RSOS model). It is straightforward to see that the limit p ^ 1 implies that a; ^ 0, since K'{p 1) 
diverges. To find out the limiting values of the weights, we use the following properties of the function E{z,x): 
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In this limit the weights of the RSOS model (unshifted as well as shifted) can be written as 
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Using the limits above and taking into account the exact forms for the composite weights (see Appendix |A]) , we are 
able to show that the diagonal weights become much larger than the off-diagonal ones. In particular, the leading x 
behavior of the weights of the diagonal plaquettes (those which are NE-SW symmetric) is given byxi-'-/2,in the limit 
X ^ 0. We explicitly checked that all the off-diagonal plaquettes have weights with a leading exponent of x strictly 
larger than 1 - r/2, showing that A is diagonal when p ^ 1. So in this limit we can drop all the off-diagonal weights 
and only diagonal weights contribute to the corner transfer matrix: 



9h9h 



(CIO) 



The dependence on w of the weights in the limit p ^ 1 follows from the form of the weights ai, etc., given in Eq. (C9 1 



combined with the form of the weights given in Eq. (C4 1. Note that the exponent of w is always integer or half-integer. 
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Inserting the above in the definition A in Eq. ( 42 1 . we get that 



(m+l)/2 

^1,1'= n [Wil2j-l,hj,l2j+l,hj+2,hj,l2j+l,l2j+2,l2j + 3)YS{l,l') (Cll) 

and 

Au-gr>^^'^'5(l,l'), (C12) 
where we have dropped the irrelevant multiplicative factor x^"'"^^ and 0(1) is given in Eq. (60 1. 



Appendix D: Connection with conformal field theory for r > 5 

In the case r > 5, we did not yet obtain explicit expressions for the functions Xm{q) with m finite. However, the 
functions Xm{q), with m odd, satisfy the following recursion relations (m is odd by definition) 

Xm{a;b,c,d,e;q) =q 2 Xm-2{a;b - 2,b - l,b,c;q) + 

q 2 Xm-2ia;b,b- l,b,c;q) + 

q " X,n-2ia;b,b+ l,b,c;q) + 

q 2 X„_2(a;&+ 2,&+ l,6,c;g) , 



(Dl) 



where we define Xm{a; b, c, d, e; g) to be zero if any of the a, b, c, d, e lies outside of the range 1, 2, . . . , r - 1. In addition, 
Xm{a]b,c,d,e]q) = if |6 - c| + 1, |c - d] + 1, \d-e\ + 1. Finally, Xi(a; 6, c, d, e; g) = 5l«-e|/4+<5a,e<5b,<i5c,. ^ if also |a _ 6| = 1 
and zero otherwise. 

Using these recursion relations, one can obtain high order expansions for Xm{q), which allows one to identify the 
conformal field theory, by direct comparison to the CFT characters. The following exact results were obtained in this 
way, and hence not proven. 

1. The case u > (domain 

For M > 0, we could identify the functions Xm{q), where {a]b,c,d^e) correspond to ground state patterns, as 
characters of the coset model ■^"(^)i''^"(^)i''^"(^)''-2 ^ -^^i^ere r - 2 = k > 3. The ground states for u > were discussed 



in section VIA The fields in this model are labeled by where t' and t correspond to the factors sm(2)i, si 

corresponds to SM(2)fc-2 and S2 to su{2)k- Because of the constraint t' + 1 + si + S2 = ^ mod 2, we can set t' = 1 (both 
t' and t can take the values 1 and 2). Finally, si = 1, 2, . . . , fc - 1 and 82 = 1, 2, . . . , A: + 1. Note that we use the height 
values to label the fields. 

To make the connection between the labels (i',t, 81,82), we will assume that m = Ap+ 3, with p an integer (the case 
m = Ap+l is very similar). Because of the relation Xm{a,', b, c, d, e; q) = Xm{r - a; r - b, r - c, r - d, r - e; q) , we only have 
to consider two cases for the labels {a;b,c,d,e). The function Imip^ao X4p+3(a;b,b - 1,6,6+ l;g) gives the character 
of the field Finally, the function limp^^o ^4p+3(a; b,b + l,b + 2,b + l;q) gives the character of the fields 

For completeness, we give the scaling dimensions of the fields in the coset theory explicitly. This formula resembles 
the formula for the scaling dimensions of the minimal models [39 . 

Finding the scaling dimensions of the fields in coset conformal field theories is typically easiest done in a Coulomb 
gas formalism HQ] . In the case at hand, the relevant Coulomb gas was studied in Ref. [TH] . Based on those results we 
find (by appropriately constraining the values of the electric and magnetic charges) 



■ (si(fe+2)-S2fe) -4 , 1 (si-S2+2t) mod 4 r . j n n 

' — — — + f - — ^"^^ lor Si + 82 mod 2 = 

+ 1 

8fe(fe+2) 8 



' — — + ^ for Si + 82 mod 2=1 



The scaling dimensions satisfy h{'i-t,k- 8i,k + 2- 82) = h{t,8i,82), reflexing the fact that the fields 'I'fe+2-s2* ^^'^ 



<j)t^,t,si a^j^g identified. 



30 



2. The case u < (domain D2) 



For M < 0, the system is described in terms of parafermions, and we find the following identification. The 
parafermion fields are labeled by two integers, where j = 0, 1, . . . , A;, and j + n = mod 2. Two fields which only 
differ in their n label by 2fc are identified, $^ e $^^2fc- addition, one has the identification = which is 

refiected in the function Xm{q) via Xm{a; b, c, d, e; q) = Xm{r - a]r - b,r - c,r - d,r - e]q). We remind the reader that 
k = r - 2, and that the parameters a, . . . , e lie in the range 1, 2, ... r - 1. 

For u < 0, the ground states are specified uniquely by (a; &, c), because d = b and e = c (the ground states for u < 
were discussed in Section VI B I. Note that c = 6 ± 1. The label a corresponds directly to the label j = a - 1. The 
magnitude of the field label n is given by c - 1, with a positive sign if c = 6 + 1, and a negative sign if c = 6 - 1. Thus, 
in general we find that liijim^oo q^"^^^^^™^^^^^Xm{a;b,c,b,c;q^^) corresponds to the character of the field ^'(^c-b){c-i)' 
where the limit is taken over the odd integers. 
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